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We present a study of the decays B° -> D ( * )0 K ( * )0 using a sample of 226 million T(45 r ) -> BB 
decays collected with the BABAR detector at the PEP-II asymmetric-energy e + e~ collider at SLAC. 
We report evidence for the decay of B° and B° mesons to the D*°K^ final state with an average 
branching fraction B(B° -» D*°K°) = (B(B° -> D*°K°) + B{B° -» D*°K°))/2 = (3.6 ± 1.2 ± 
0.3) x 10~ 5 . Similarly, we measure B(B° -> D°K°) = (B{B^_ -> D°J^) + 2?(B° -> D°K°))/2 = 
(5.3±0.7±0.3) x 10~ 5 for the D°K° S final state. We measure B(B° -> D°K*°) = (4.0±0.7±0.3) x 10~ 5 
and set a 90% confidence level upper limit B(B° — > D /^* ) < 1.1 x 10 -5 . We determine the upper 
limit for the decay amplitude ratio \A(B° -> _D°_K r *°)/yl(B -» D°F*°)| to be less than 0.4 at the 
90% confidence level. 



PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 



With the discovery of CP violation in the decays of 
neutral B mesons [lj and the precise measurement |2j of 
the angle j3 of the Cabibbo-Kobayashi-Maskawa (CKM) 
Unitarity Triangle (3j , the experimental focus has shifted 
towards over-constraining the unitarity triangle through 
precise measurements of \V u i>\ and the angles a and 7. 
The angle 7 is axg(-V; b V ud /V* b V cd ) and V l} are CKM 
matrix elements. Several methods have been suggested 
and explored to measure 7 with small uncertainties |^|, 
but they all require large samples of B mesons not yet 
available. The decay modes B° — > D^°K° offer a new 
approach for the determination of sin(2/3 + 7) from the 
measurement of time-dependent CP asymmetries in these 
decays [B| . The CP asymmetry appears as a result of the 
interference between two diagrams leading to the same 
final state D^K® (Figure [IJ. A B° meson can either 
decay via afe^c quark transition to the D^*>°K° (K° — ► 
Kg) final state, or oscillate into a B° which then decays 
via a b -> u transition to the D^K {K° -> K°) final 
state |6(. The J5°_B° oscillation provides the weak phase 
2/3 and the relative weak phase between the two decay 
diagrams is 7. 



D (*)o 




K (*)o 



FIG. 1: The decay diagrams for the b — *■ c transition B — > 
D ( * )0 K° and the b -> u transition B° -> D ( * )0 K°. 

The sensitivity of this method [5| depends on the rates 
for these decays and the ratio of the decay amplitudes. 
The branching fractions B(B° -* Z)(*) J?(*) ) can be esti- 
mated from the measured color-suppressed decays B° — ► 
£ > (*)o 7r o j^ t0 be approximately B(B° -> £)(*) j?(*) ) w 
sin 2 6> c B(B° -> DM°7r ) - e>(10~ 5 ), where C is the 
Cabibbo angle and sin^c = 0.22. The Belle Collabora- 



tion has observed the B Q —> D°K l -*' )t) decays with branch- 
ing fractions consistent with this naive expectation y|. 
The time-dependent CP asymmetries in B° — > D^*'°K° 

decays are proportional to r B •sin(2/3 + 7±o')/(H-rg ), 



where 



„(*) - 



\A{B° -> D^°K )/A(B° -» Df*) ^ )! 
and <5 is a relative strong phase which depends on the 
specific final state. Higher values of r B lead to larger 
interference between the b — > c and b — > u processes and 
thus increased sensitivity to the angle 7. In the Standard 



Model r 



(*) 



/ ' \V u bV* s \/\VcbV* s \, where the factor / ac- 



counts for the difference in the strong interaction dynam- 
ics between the b — > c and 6 — » u processes. There are no 
theoretical calculations or experimental constraints on /. 

In B° -v D(*)°F (JT° ->■ AT°) decays the strangeness 
content of the if is hidden and one cannot distinguish 
between B° -> D^ a K Q and B° -> D^ a K a . Therefore 
a direct determination of r B from the measured rates is 
not feasible. In the remainder of this paper we refer to 
these decays as B° -> D^K . Insight into the B de- 
cay dynamics affecting r B can be gained by measuring a 
similar amplitude ratio r B = \A(B° -> D°K*°)/A{B° -> 
£)°A"*°)| using the self-tagging decay if* — * A~7r + . The 
B° -» £)°F*° and _B° -+ ZJ°F*° decays are distinguished 
by the correlation between the charges of the kaons pro- 
duced in the decays of the neutral D and the K*°. In the 
former decay the two kaons in the final state must have 
the same charge, while in the latter they are oppositely 
charged. This charge correlation in the final state is di- 
luted by the presence of the doubly-Cabibbo-suppressed 
decays D° — ► K + tt~ ,K + ir~ir°, and K + ir~Tr + ir~ . The 
ratio Tb is related to the experimental observables IZi 
defined as 



Hi 



T(B° -* (K+X7) D K*°) 
r(so - (K-Xf) D K*°) 

?b + r D, + % r ' B r Di cos(7 + Si), 



(1) 



where 



Xf = 

r D, = 

6, = 



7T ,7T 7T ,7T 7T 7T + , 

\A(D°^K+Xr)\ 
\A(D°^K-Xf)\' 

5b + 5 Di , 



(2) 
(3) 

(4) 



and 8b and Sdi are strong phase differences between the 
two B and Di decay amplitudes, respectively. The values 
of rDi have been measured to be Td^k-k = 0.060 ±0.002, 
r D ^K^» = 0.066 ± 0.010, and t D ->k-k-k-k = 0.065 ± 
0.010 0. 

We present herein measurements of the branching frac- 
tions B(B° -y L> A°) and B(£° -> D°K*°), evidence for 
the decay B° -> D*°K°, a 90% confidence level (C.L.) 
upper limit for the branching fraction of the b — > u tran- 
sition B° — > D°K*°, and a limit for the ratio fs- 

These results are based on a sample of 226 million 
T(45) -> SB decays collected with the BABAR detector 
between 1999 and 2004 at the PEP-II asymmetric-energy 
e + e~ collider operating at the T(4S) resonance. The 
properties of the continuum e + e~ —* qq (q — u, d, s, c) 
background events are studied with a data sample of 
11.9 fb _1 recorded at an energy 40 MeV below the 
T(4S) resonance. The BABAR detector has been de- 
scribed in detail elsewhere [lfj. Detector components 
relevant for this analysis are summarized here. Trajec- 
tories of charged particles are measured in a spectrome- 
ter consisting of a five-layer silicon vertex tracker (SVT) 
and a 40-layer drift chamber (DCH) operating in a 1.5 T 
axial magnetic field. Charged particles are identified as 
pions or kaons using information from a detector of inter- 
nally reflected Cherenkov light, as well as measurements 
of energy loss from ionization (dE/dx) in the SVT and 
the DCH. Photons are detected using an electromagnetic 
calorimeter composed of 6580 thallium-doped Csl crys- 
tals. We use a Monte Carlo simulation of the BABAR 
detector based on GEANT4 [llj to validate the analy- 
sis procedure and to study the backgrounds. Simulated 
events are generated with the EvtGen |l2( event genera- 
tor. 

We reconstruct the decays B° -> D°K°, D*°K°, 
D°K*°, and D°K*° in the decay chains: D^ -> L>°7r°; 
D" -> K~ir+, K-tt+tt , and K-ir+ir-ir+; K° -> K° -> 
7r + 7r~; K*° — » K~ir + ; and ir° — > 77. For each B de- 
cay channel the optimal selection criteria are determined 
by maximizing the ratio Ns/-\/Ns + Nb , where Ns and 
Nb are, respectively, the expected signal and background 
yields estimated from samples of simulated events. A 
large sample of the more abundant B + — ► D°ir + de- 
cays, in which the D° decays to the K + ir~ , K + tt~it , or 
K + ir~ir + Tr~ final states, is used as a calibration sample 
to measure efficiencies and experimental resolutions for 
the selection variables. 

Well reconstructed charged tracks are used to recon- 
struct D° and K*° candidates. The K^ candidates must 



satisfy a set of kaon identification criteria. 

These identification criteria have an average efficiency 
of about 90%, while the probability of a pion to be 
misidentified as a kaon varies between a few percent and 
15%. Photons are reconstructed from energy deposition 
clusters in the electromagnetic calorimeter consistent 
with photon showers, and are required to have an energy 
greater than 30 MeV. We select tt° candidates from pairs 
of photon candidates by requiring their invariant mass to 
be in the interval 115 MeV/c 2 < 01(77) < 150 MeV/c 2 . 

The K° candidates are selected from pairs of op- 
positely charged tracks with invariant mass within 
7 MeV/c 2 (~ 2cr) of the nominal K° mass. The dis- 
placement of the K® decay vertex from the interaction 
point, in the plane perpendicular to the beam axis, di- 
vided by its estimated uncertainty must be greater than 
2. The K*° candidates are selected from pairs of oppo- 
sitely charged K + and 7r _ tracks, with invariant mass 
within 50 MeV/c 2 of the nominal K*° mass. The polar- 
ization of the K*° in the B° decay is used to reject back- 
grounds by requiring |cos^^| > 0.4, where the helicity 
angle Oh is defined as the angle between the direction of 
the K*° in the B° meson rest frame and the direction of 
its daughter K+ in the K*° rest frame. For B° -» D°K*° 
and B° — » D°K*° signal candidates, 9h follows a cos 2 9h 
distribution, while the combinatorial background is dis- 
tributed uniformly. 

We reconstruct D° candidates in the K~n + and 
A'~7T + 7r~7r + decay modes by combining charged tracks, 
retaining combinations with an invariant mass within 2a 
of the nominal D° mass m^o. In the D° — ► K~ir + Tr° se- 
lection, the 7T° candidates are required to have a center- 
of-mass (CM) momentum p* greater than 400 MeV/c. 
For each A~7r + 7r° combination, we use the kinematics 
of the decay products and the known properties of the 
Dalitz plot for this decay [P3j to compute the square of 
the decay amplitude A 2 . We select combinations with 
A 2 greater than 5% of its maximum value. This require- 
ment selects mostly the K~ p + region of the Dalitz plot. 
It rejects 62% of the combinatorial background, while 
keeping 76% of D° — ► K~tt + tt° signal, as measured with 
the B + — ► D°7r + sample. Combinations with invariant 
mass within 25 MeV/c 2 (2. 5c) of m D o are retained. 

The D*° candidates are selected from combinations 
of a D° and a tt° with p* n0 > 70 MeV/c. After kinc- 
matically constraining D° and ir° candidates to their 
nominal masses, we select the candidates with a mass 
difference Am = |m(L>*°) - m(D°) - 142.2 MeV/c 2 | < 
3.3 MeV/c 2 (3c). 

Two standard kinematic variables are used to select 
B° candidates: the energy-substituted mass mgsc 2 = 



v^ 



: 2 pr • Pb) 2 /E 2 - — c 2 p 2 B and the energy differ- 
ence AE = E B — |-/s, where the asterisk denotes the 
CM frame, s is the square of the total energy in the CM 
frame, p and E are, respectively, three-momentum and 



energy, and the subscripts T and B refer to T(4S) and 
B . In calculating p# and E* B we constrain the mass of 
the D(*>° and K® candidates to their respective nominal 
values. For signal events, toes is centered around the B° 
mass with a resolution of about 2.6 MeV/c 2 , dominated 
by knowledge of the e + and e~ beam energies. In simu- 
lated events the AE resolution is found to be w 13 MeV 
for all B° decay modes considered in this analysis. The 
B° candidates are required to have Toes > 5.2 GeV/c 2 
and \AE\ < 100 MeV. 

We use two variables to reject most of the remaining 
background, which is dominated by continuum events: 
a Fisher discriminant [14j based on the energy flow in 
the event and the polar angle 9 B of the B° candidate 
in the CM frame. For correctly reconstructed B candi- 
dates cos 9 B follows a 1 — cos 2 9 B distribution, whereas it 
is uniformly distributed for continuum events and com- 
binatorial background. We require |cos6*J| < 0.75 for 
B° -» D Q K*°, and |cos6^| < 0.85 for all other decay 
modes. The Fisher discriminant T is defined as a linear 
combination of | cos 9^ B | and two energy- flow moments 
C° and C 2 . The variable 9^ „ is the angle in the CM 



frame between the thrust axis |15j of the decay products 
of the B° and the thrust axis of all charged and neutral 
particles in the event excluding the ones that form the 
B°. The energy-flow moments £° and C 2 are defined as 
C l = J2jP*j cos * 9j where p* is the CM momentum and 
9j is the angle between the direction of particle j with 
respect to the thrust axis of the B° candidate, and the 
sum is over all particles in the event (excluding those 
that form the B°). The requirement on T varies for 
each decay channel because of different levels of expected 
background. In the D^°K° and D°K*° final states our 
requirement has an efficiency of about 80% for the signal 
while rejecting approximately 85% of the background; in 
the B° — > D°K*° mode a tighter requirement rejects 95% 
of the background and has a signal efficiency of 55%. 

In the D*°K® final state, approximately 5% of the 
events that satisfy all selection criteria contain more than 
one B° candidate. We retain the candidate with the 
smallest % 2 computed from the measured value of m(D a ) 
and m(D*°) — to(Z?°), their nominal values, and their 
resolutions in data. In the D°K° S , D°K* a , and D°K*° 
final states we retain all selected B° candidates since the 
fraction of events with two or more candidates is negligi- 
ble (< 1%). 

The selected B° ->■ D^K^ candidates include 
small contributions from numbers of B decays to simi- 
lar final states which are misreconstructed as signal can- 
didates. We have studied these backgrounds with large 
samples of simulated events, corresponding to between 
100 and 1000 times the size of our data sample, for the 
following categories of decays: (1) B° —» D°p°,p° — > 
7r + 7r~ decays, where one of the two pions is misiden- 
tified as a charged kaon; (2) B° — » D + ir~ decays fol- 
lowed by Cabibbo-suppressed decays D + — > K^*' Q K + , 



and B° -* D+R- followed by D+ -> K^°ir + , recon- 
structed in the D a (K-ir + )K ( -*^ final states; (3) charm- 
less B° — ► K~Tr + K®(nir) where the K~ and 7r + are 
wrongly combined to form a D° — > K~ir + candidate; 
(4) B a -> D*°if(*)°, D*° -* _D°7 candidates, where a 
low energy photon is not reconstructed; (5) the decays 
B~ -> D*°K-,D*° -> D°tt°/j, B~ -> D°K*-,K*- -> 
K-tt°,K° s tt-, and B° -> D*-R+,D*- -» D°ir- , where 
a low-energy 7r°, 7r _ , or photon is replaced by a random 
low-momentum charged particle. The contribution of 
category (1) is found to be less than 0.01 events and hence 
is neglected. The contribution of category (2) is also neg- 
ligible in all modes, except for B -> D°K°, D° -> K~tt + . 
We eliminate 87% of these background events by requir- 
ing the invariant masses m(K®K + ) and m(K^TT + ) to be 
more than 20 MeV/c 2 away from the nominal D + mass. 
The toes spectrum of the remaining background events 
in this category, and in categories (3)-(5), show a broad 
enhancement near the B mass. However, due to the D° 
mass constraint, B° candidates with a misreconstructed 
D° do not peak, unlike the signal, in the AE distribution 
at zero. In the decay B° — > D°K*°, the charge correla- 
tion used in the selection removes all contributions from 
known B decays included in simulated events. 

The signal yield for each B° decay mode is deter- 
mined with a two-dimensional extended unbinned maxi- 
mum likelihood fit to the Toes and AE distributions, sep- 
arately for each D° decay mode. The probability density 
function (PDF) is a sum of three components: a signal 
component (5(toes) x Q{AE), a background component 
G(m-Es) x 'Pi(AE), accounting for other B decays mis- 
reconstructed as signal, and a combinatorial background 
component T(toes) X Vi{AE). Here, <7(toes) is a Gaus- 
sian describing the toes distribution of signal and mis- 
reconstructed B decays; Q(AE) is a Gaussian describing 
the signal AE distribution; T'i(AE) are first-order poly- 
nomials describing the AE distributions of background 
events. The toes distribution of the combinatorial back- 
ground is parameterized by a threshold function T(toes) 
defined as T(toes) ~ "Iesa/1 — x 2 exp{— £(1 — x 2 )} [l6|. 
where x = 2toesc 2 / \fs and £ is a shape parameter. The 
mean and the resolution of Q (toes) and Q(AE) are fixed 
to values measured in the B + — ► D°n + calibration sam- 
ple. 

The measured signal yields are summarized in Table [I] 
The AE distributions of candidates with |toes — 5280| < 
8 MeV/c 2 for the sums of the reconstructed D° decay 
modes are illustrated in Figure |2J The signal significance 
S is computed as S = v /2(ln£(iV s ) - ln£(N s = 0)), 
where C(Ns) is the maximum likelihood of the nominal 
fit, and C(Ns = 0) is the value obtained after repeating 
the fit with the signal yield N$ constrained to be zero. 

The branching fraction B for each B° decay mode is 
the weighted average of the branching fractions Bj in 
each D° channel Dj = {K~ir + , K~ir + Tr~n + , K~tt + tt }, 



TABLE I: Signal yield Ns, signal significance 5, effective 
signal efficiency e e g, and the measured branching fraction B 
for the B° -► D^°K°, B° -> D°K*°, and B° -> _D°F*° 
decays. The efficiency e c fr is defined as £^\ e$ x Bi, where the 
sum is over the D decay modes, e% are the signal reconstruc- 
tion efficiencies, and Bi are the corresponding intermediate 
branching fractions for D*°, D°, K*°, and K° decays to final 
states reconstructed in this analysis. 



B Mode 



Ns 



S e c « [%] B [10- 



B° -* D°K° 104 ± 14 9.2a 
D*°K° 17.1 ±5.2 4.3a 



i) 



B 

B° -► B°r° 

5° -► ^°f*° 



77 ±12 7.9cr 



-3.6 



+6.8 



0.82 
0.17 
0.84 
0.47 



5.3 ±0.7 ±0.3 
3.6 ±1.2 ±0.3 
4.0 ±0.7 ±0.3 
0.0 ±0.5 ±0.3 
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FIG. 2: Distribution of AE for a) £" 

D* ^ , c) B° -> D°K* , and d) B° -*■ _D°A"* U candidates 
with |m ES - 5280 MeV/c 2 | < 8MeV/c 2 . The points are the 
data, the solid curve is the projection of the likelihood fit, and 
the dashed curve represents the background component. 



computed as 



B, 



N Sl 



2 x N BE x B(T(45) -» S°S°) x B Dj x6 k xe. 



■(5) 



where Ns- is the signal yield from the likelihood fit, 
N BB - is the total number of T(4S I ) — ► SB events, Bd, 
is the branching fraction B(D° -> D,-) in B° -> D Q K^° 



in B° -> D* ^ , 
-> i4T + 7r _ ) branch- 
► D K*°,D°K* ), 



and £(L>*° -» L>°7r°) x B(D° -> D 6 
B K is the K° -> K° s -> 7 r+7r^( J r£*° 
ing fraction in 5° -> D^°K (B° - 
and £j is the signal reconstruction efficiency. We as- 
sume B(T(4S) -* -B°B°) = 0.5. The systematic un- 
certainties for the branching fractions include contribu- 
tions from estimated misreconstructed B background (1- 
13%) |l7f . variation of parameters kept fixed in the 
likelihood fit (2-8%), D^° branching fraction (2.4- 



6.9%), 7T° reconstruction efficiency (3%), photon re- 



construction efficiency (1.8%), charged-track reconstruc- 
tion efficiency (0.8% per track), simulation statistics (1- 
4%), efficiency correction factors (1-4%), kaon identifi- 
cation efficiency (2% per kaon), K® reconstruction ef- 
ficiency (1.6%), and the number of BB events (1.1%). 
The efficiency correction factors are obtained by com- 
paring data with MC simulation in the B + — ► D°ir + 
control sample. The largest contributions to the uncer- 
tainties in these factors are from selection requirements 
for the 7T° momentum p* and the amplitude |„4| 2 in the 
decay and the Fisher discriminant J- ' . 



We measure 



B(B° -> D°K°) 
B(B° -> D*°K°) 
B(B° -► D°K*°) 
B(B° -» D Q K*°) 



(5.3 ±0.7 ±0.3) x 10~ 5 
(3.6 ±1.2 ±0.3) x 10~ 5 
(4.0 ±0.7 ±0.3) x 10~ 5 
(0.0 ±0.5 ±0.3) x 10~ 5 



where the uncertainties are, respectively, statistical and 
systematic. For the decay B° — ► D°K*° we use the 
Bayesian method to compute the upper limit Nul on 
the observed number of events. The value of Nul at 90% 
C.L. is defined as f* UL C(N) dN = 0.9, where C(N) is 
the maximum likelihood function from the fit to the tties 
and AE distributions. We assume a flat prior probabil- 
ity density function for B > 0. We account for system- 
atic uncertainties by numerically convolving C(N) with a 
Gaussian distribution with a width determined by the rel- 
ative systematic uncertainty multiplied by the measured 
signal yield. We obtain B(B° -» D°K*°) < 1.1 x 10~ 5 
at 90% C.L. 

We compute an upper limit on the ratio rg by mea- 
suring the ratio TZi in each D° decay mode. We use 
the expression TZi = (e D kI e d k) ' (^dk/^dk) to ^" 
tain the PDF for TZi from the unbinned maximum likeli- 
hood fit described earlier. In this expression £j) k ( £ d k) 
and Nqk {N D ^ ) are, respectively, the reconstruction 
efficiency and fitted _yield of the B° -> D°K*°,D° -> 
K~X? (B° -» D°K*°,D Q -* K-X+) decay modes. 



The uncertainties on e/jj^ , 



e D Ri ano - ^d k are use d 



to obtain the posterior PDF CiTZi) for each IZi. We as- 
sume a Gaussian PDF for ro i . We compute the PDF for 
fs by convolving C(lZi) and rn i according to equation 
0. We obtain the limit f B < 0.40 at 90% C.L. with a 
Bayesian method using uniform priors for TZi > and 
by taking into account the full range 0°-180° for 7 and 
Si- The present signal yields combined with this limit 
on ?b suggest that a substantially larger data sample is 
needed for a competitive time-dependent measurement 
of sin(2/3 + 7) in B° -+ D^K decays. 

In summary, we have presented measurements of the 
branching fractions for the decays B° — ► D°K° and 
B° -> D°K*°, evidence for the decay B° -> D*°K Q , and 
an upper limit for the ratio tb ■ Our results are in agree- 
ment with previous measurements of these modes 18]. 
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